Condensation of actin filaments pushing against a 
barrier 



Kostas Tsekouras^'^, David Lacoste^, Kirone Mallick^ and 
Jean-FranQois Joanny^ 

1 Laboratoire de Physico-Chimie Theorique - UMR CNRS Gulliver 7083, ESPCI, 10 

rue Vauquelin, F-75231 Paris, France 

^ Physico-Chimie UMR 168, Institut Curie, Paris, France 

^ Service de Physique Theorique, Commissariat a I'Energie Atomique- Saclay, Gif, 
Prance 

E-mail: tamytes8a@gmail . com 

Abstract. We develop a model to describe the force generated by the polymerization 
of an array of parallel biofilaments. The filaments are assumed to be coupled only 
through mechanical contact with a movable barrier. We calculate the filament density 
distribution and the force-velocity relation with a mean-field approach combined with 
simulations. We identify two regimes: a non-condensed regime at low force in which 
filaments are spread out spatially, and a condensed regime at high force in which 
filaments accumulate near the barrier. We confirm a result previously known from 
other related studies, namely that the stall force is equal to N times the stall force of 
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the velocity is practically zero at forces significantly lower than the stall force. 
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1. Introduction 

Actin filaments and microtubules are key components of tlie cytoskeleton of eukaryotic 
cells. Both play an essential role for cell motility and form the core components of 
various structures such as lamellipodia or filopodia. They are active elements which 
exhibit a rich dynamic behavior. For instance, actin filaments treadmill in a process 
where monomers are depolymerized from one end of the filament while other monomers 
are repolymerized at the other end. Actin polymerization is highly regulated in the 
cell, through many actin binding proteins. Some of these proteins accelerate actin 
polymerization, while others crosslink filaments or create new branches from existing 
filaments. All these proteins ultimately control the force that a cell is able to produce 
pp. 

Given the complexity of actin polymerization, many studies have focused on its 
basic structural element, namely the filament itself. For instance, a lower bound for 
the polymerization force generated by a single actin filament has been deduced from 
the buckling of a filament which was held at one end by a formin domain and at the 
other end by a myosin motor P]. Other studies focused on the dynamics of single 
filaments through depolymerization experiments [3]. In order to understand the rich 
dynamical behavior of single filaments like actin or microtubules and the force they 
can generate, discrete stochastic models have been developed which incorporate at the 
molecular level the coupling of hydrolysis and polymerization [U |5l [6l [71 [H HI [TOl [11] . 
The filament dynamics and the force generation are two related aspects: Hydrolysis is 
not only relevant for understanding the single filament dynamics but also for the force 
generation, since the force generated by a filament is typically lowered by hydrolysis [0] . 

Ensembles of parallel interacting filaments are able to generate larger forces than 
single filaments as in cellular structures called filopodia [12]. General thermodynamic 
principles controlling the force produced by the polymerization of growing filaments 
pushing against a movable barrier were put forward many years ago by Hill et al. 
[13]. For many years however, it was unclear how to extend these results in order 
to understand theoretically the effect of interactions or collective effects in the process 
of force generation. Progress in this direction was made through the introduction of 
stochastic models for ensembles of parallel microtubules [HI [IHl [16], and through the 
development of simulations for actin filaments in parallel geometry [17] or in networks 
[IS] . In these works, the brownian ratchet model [T^] was used at the single filament 
level, while some specific rule was assumed concerning the way the load is shared between 
the filaments. In the absence of hydrolysis and lateral interactions between the filaments, 
the stall force of an ensemble of parallel filaments should be A^ times the stall force 
of a single filament, as confirmed either by a detailed balance argument valid only near 
stalling [15] or more recently by a more general analysis based on a decomposition into 
cycles [2b]. 

In this work, we propose a new theoretical framework for this problem. One 
novel aspect as compared to previous work [15j comes from the fact that we model 
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the dynamics of an ensemble of parallel non-interacting filaments at an arbitrary value 
of the force, rather than just predict the value of the stall force. Another important 
difference between this model and previous work, is that our model allows an arbitrary 
number of filaments in contact with the moving wall, which allows the possibility of a 
condensation transition for the number of filaments at the wall. 

This paper is organized as follows: we first present the model, secondly the mean- 
field approach for the general case of an arbitrary A^, then the simulations, and a 
theoretical analysis of the approach to stalling. We end with a discussion of various 
related experiments in this field, in which forces generated by a few actin filaments have 
been measured [22] . 

2. Model 

We consider two rigid fiat surfaces: one fixed where filaments are nucleated (nucleating 
wall) and one movable (barrier) whose position is defined to be the position of the 
filament (s) furthest away from the nucleating wall (thus there is always at least one 
filament in contact with the barrier). In the cellular environment, this "barrier" is 
often a membrane against which filaments exert mechanical forces. We do not model 
the internal structure of the filaments, and in particular we do not account for ATP 
hydrolysis. After nucleation, the filaments grow or shrink by exchanging monomers 
with the surrounding pool of monomers, which acts as a reservoir. The filaments are 
coupled only through mechanical contact with the barrier. In some previous models 
pH], a staggered distribution of initial filaments was assumed so that there would be 
only a single filament in contact at a time. Here we do not make such an assumption, 
on the contrary monomers inside different filament are precisely lined up. As a result, 
the number of filaments at contact is an arbitrary strictly positive integer. 

It follows that we can separate the filaments in two populations, the free filaments 
which are not in contact with the barrier, and the filaments in contact. Only the 
filaments in contact feel the force exerted by the barrier on them, and as a result this 
changes their polymerization rates as compared with free filaments. We assume that a 
monomer can be added to any free filament with rate Uq or removed with rate Wq, as 
shown in Fi^ Similarly, a monomer can be added to a filament in contact with rate 
U{F), and removed with a rate W{F) (or Wq as explained below). The values of the 
rates which we have used correspond to an actin barbed end and are given in table [T] 
We also assume that the barrier exerts a constant force F on the filaments in contact, 
this force is defined to be positive when the filaments are compressed. 

We need now to specify more precisely how the force exerted by the barrier is 
shared by the filaments in contact. When a monomer is added to a filament in contact, 
the barrier moves by one unit, but only the filament on which the monomer has been 
added does work; we therefore treat all the other filaments as free during that step. 
Similarly, during depolymerization, filaments depolymerize from the barrier with the free 
depolymerization rate Wq as long as there is at least one other filament in contact with 
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the barrier, since in this case the depolymerizing filaments do not produce work. The 
depolymerization occurs with a rate W only when there is a single filament in contact 
with the barrier. In this case the filament produces work, since its depolymerization 
leads to the motion of the barrier. 

For a filament which has exchanged work with the barrier through addition or loss 
of monomers, we use a form of local detailed balance which reads : 

^ = ^e-/. (1) 

W Wo ^ ^ 

This relation is obeyed by the following parametrization of the rates [T^ |21]: 

U = U^e-f^ and W = Woc^^^'^K (2) 

where 7 is the "load factor" and / is the non-dimensional force / = Fd/ksT, where d 
is the monomer length. Note that 7 itself could be a function of the force, however in 
the following we assume that it is just a constant. More elaborate treatments of the 
load dependence of the transition rates can be found for instance in Ref. |25j . 

An essential feature of this model is that although multiple filaments interact with 
the barrier, when a monomer is added to one of the filaments in contact, it must do 
work against the entire load. In the classification of [26], this corresponds to a scenario 
with "no load sharing". If the force could be shared by more than one filament or if 
the monomers in different filaments were not precisely lined up, the above discussion 
would still apply: in this case a single filament would carry a fraction of the load at a 
time, and for that filament a similar local detailed balance would hold. In this case, 
although the stall force would be the same as in the "no-load sharing" scenario, the 
form of the force- velocity curve would be affected. Such models have been considered in 
Refs. [15], [ini EHl [20], but for simplicity, in the present paper, we focus on the "no load 
sharing" model. 

3. Theory 

In the particular case that there are only two filaments (A^ = 2), the master equation 
can be solved exactly in terms of the probability that there is a given gap at a given 
time between the two filaments, as shown Appendix A. Unfortunately, this approach 
is limited to the N = 2 case, because only in that case there is a single gap between 
the filaments. For N > 2, there are many gaps, so in general such an approach quickly 
becomes as complicated as the one based on the filaments themselves. So instead of 
looking for an exact solution, we provide in the section below, an approximate but 
accurate mean- field solution for the general case N > 2. 

3.1. An ensemble of N filaments with N > 2 

We recall that the position of the moving barrier coincides with the position of the 
longest filament, and we define A^, as the number of filament ends, which are present at 
a distance i from the moving barrier. We take the convention that i = corresponds 
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to the barrier itself. Since each filament has only one active end and the total number 
of filaments is fixed to be N, we have the condition that X]i=o ^* ~ ^ft^r a careful 
account of all the possible events that can occur on any filament in a small time interval, 
we obtain the following master equations: 

"^^^ {Wo+UNo)N,^,+{Uo+W5mo=i)N,+,-{Wo+Uo+UNo+WSn,=i)N„{3) 



dt 
dNi 

~dr 

dNo 



(?7o+Ty(5;vo=i)A^2-(f/o+W^o+W^^7Vo=i+?7Aro)Ari+[l^o(l-57Vo=i)+f/(A^o-l)]iVo,(4) 
{Uo + W6r,,=i)N, - [U{No - 1) + 1^0(1 - Sno=i)]No, (5) 



dt 

where 6no=i represents the probability that there is only a single filament in contact. 

In deriving these equations, we have for instance implicitly replaced the joint 
probability to have Ni filament ends at position i and to have only one filament at 
contact at time t, namely P{Ni{t) = Ni,No{t) = 1) by the product of P{N,{t) = Ni) 
and P{NQ{t) = 1). In other words, a mean-field approximation has already been used. 
A further consequence of this mean-field approximation is that in these equations, 6no=i 
can be replaced by its time-averaged value, which we call q: 

q = {Sno=i)- (6) 

The quantity g is a central feature of our model for N > 2. All subsequent results and 
calculations appearing in this paper follow from this mean-field approximation. 

At steady state, the l.h.s. of Eq. [3] is zero. The r.h.s. leads to a recursion valid for 
i > 2, which can be solved after a few lines of calculations. The solution is 

N, = N2exp{-{t-2)/l), (7) 

where / is the correlation length (expressed in number of subunits) given by 

Uo + Wq " 
Wo + UNo 

The other two equations Eqs. |4]j5l together with the normalization condition fix A''2, A^i 



In 



(8) 



and Nq. We find that the average number of filaments in contact with the wall A"o is: 

= (^0 + ^1- ^0)^ .9) 

Uo + U{N-l) + {W-Wo)q ^ ' 

When N = 2, this mean-field solution agrees with the exact solution derived in 
appendix A with the additional condition that 7 = 1, in which case the on-rate carries 
all the force dependence. For an arbitrary value of 7, the mean-field solution does not 
agree with the exact result obtained for N = 2. This is expected since the mean-field 
approximation should work well only in the limit of large N. 

The average velocity of the moving barrier is 

V = d{UNo - Wq), (10) 

where the first term within the parenthesis is the contribution of the filaments in contact 
polymerizing with rate U and the second term is the contribution from depolymerizing 
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events of a single filament in contact. We have not found a way to solve in general 
the self-consistent equation satisfied by q, namely Eq. [6} except near stalling conditions 
as explained in the next section. For this reason, we have calculated numerically q 
from simulations, and derived predictions from the mean-field theory assuming that q 
is known. For instance, using Eqs. |9]and [10} one obtains the average velocity. 



4. Results 



4.1. Numerical validation of the mean-field approach 

We have tested the validity of the mean-field approach using numerical simulations. We 
used the classical Gillespie algorithm [27] incorporating the Mersenne Twister random 
number generator. Runs were executed for up to 5000. Up to 200 trial runs were used 
to derive averages and distributions. We validated the simulation results by comparing 
them with the particular cases = 1 and for N = 2 for which an exact solution is 
known (it is given in [^ for = 1 and in the previous section for N = 2). 

By evaluating the parameter q from the simulations, we obtained a very good 
agreement between the theoretical approach based on the use of mean-field and the 
simulations for the determination of the force velocity curve (shown in Fig. |2| bottom) 
and for the number of filaments A^^o in contact with the barrier (shown in Fig. [2| top). We 
find that the values of Ni as determined by theory does not deviate from the simulation 
value by more than one. 



4-2. Condensation transition as function of the applied force 

At low forces, the barrier velocity is close to its maximum value given by the free 
polymerization velocity. In this case, only one or a small number of filaments are 
in contact, therefore g ^ 1, which corresponds to a non-condensed or single filament 
regime. The steady state density profile of the filaments is broad as shown in Fig. [2] 
(bottom, left inset) and the corresponding correlation length is large. With the 
parameters values corresponding to this figure, we have / ~ 151nm. 

Inversely, at high forces, the filaments accumulate at the barrier. As a result q — 0, 
the density profile is an exponential as shown in Fig. |2[ (bottom, right inset) with a 
very short correlation length of the order of a monomer size. With the parameters 
values corresponding to this figure, we have / ~ 4.1nm. Since in this case, the number 
of filaments in contact, A"o is a finite fraction of A^, we call this regime the condensed 
regime. In this high force regime (typically near the stall force F = Fstaii)-, the following 
condition is obeyed: NU <^ Uq. Since we also have q — 0, Eq. [9] simplifies to 

No=(l- ^) N. (11) 



This equation can be used to predict the finite fraction of filaments in contact in the 
condensed regime. This condensed regime corresponds to the plateau in the curve of 



A'o vs. F which is shown in Fig. M (top inset). In the conditions of this figure, Eq. 11 
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predicts a plateau for A^o ~ N/2 = 50 which is indeed observed, and as expected the 
plateau in Nq (Fig. |2| top) occurs at the same force at which the velocity approaches 
zero (Fig.|2| bottom). 



4-3. Theoretical stall force 



Let us first discuss here the theoretical expression of the stall force and then in the 
next section the practical way this limit is approached. The stall force is defined as the 



value of the force applied on the barrier for which the velocity given by Eq. 10 vanishes. 
For iV = 1, the stall force is F^}}„ = kBT\n{Uo/Wo)/d. For = 2, using the results 

.(2) 



stall 

obtained in Appendix A for A'o and q, we find that the stall force F^^^^, is exactly twice 



the stall force of a single filament, F^^^^, 



7(2) 

stall 



2F, 



(1) 



stall 



d 



Wo 



(12) 



In the general case of an arbitrary number of filaments A^, we expect that stall force 

mm- 



F^tl should be 



stall 



d Wo 



(13) 



This result can be derived from the following argument: near stalling conditions. 



the average density of filaments at contact Nq/N, can be obtained from Eq. 11 above. 
This average density of filaments can be used as an approximation of the probability to 
have one filament in contact when Nq/N ^ 1. Since q is the probability that there is 
a single filament in contact (in other words, there is one filament among A^ in contact 
and the remaining A^ — 1 are free) , it follows that 

1 / A^ I A^ 



(14) 



which leads using Eq. [IT] to 

Wo 
Uo 



N 1 



Wo 
Uo 



N-l 



No 



Wo 

Uo 



N-l 



(15) 



We call this the binomial form for q. We note that Eq. 14 also means that 



g~ Aroexp(-Aro), (16) 
which corresponds to a Poisson statistics for the distribution of the number of filaments 



at contact. Now inserting the final expression for q of Eq. 15 into the stalling condition 



namely the vanishing of the velocity given by Eq. 10, one obtains the theoretical stall 



force given in Eq. 13 



The theoretical expression of the stall force given by Eq. 13 has been also obtained 
in a recent study devoted to the stall force of a bundle of filaments [20]. This study is 
based on the model introduced in Refs. [SI US] which the authors modified to include 
lateral interactions between the filaments of the bundle. Using a theoretical argument 
based on the identification of relevant polymerization cycles, the authors of Ref. 
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confirm the expression of the stall force obtained before in [15], which is also our Eq. 13 
More importantly, they show with this method that this expression has a universal 
character for models of this kind, hence in particular the independence of the stall force 
with respect to the load distribution factor 7. They also obtained force velocity curves 
for various values of the lateral interaction and staggering distance, which -as we have 
checked- agree with the numerical results obtained in this paper, when there is no lateral 
interaction and when the shifts are zero. 

In Fig. [3} the value of q determined from the simulations is compared with 



theoretical expression given by Eq. 14 or Eq. 16 (both expressions give similar results). 
We note that the deviation between the simulation points and the theory increases as 
the force is lowered, this is due to the mean-field nature of the theory which becomes 
invalid when the force is small since then the fluctuations are large. For completeness, 
we also show in Fig. |4] the probability density function of the number of filaments at 
contact for various forces. 

4.4- The approach to stalling 

Let us now discuss more precisely how the velocity approaches zero. We find that 
in our simulations, for N larger than about 10, the velocity approaches zero at forces 
significantly lower than the stall force as shown in Fig|2] (bottom). We note that a similar 
effect has been obtained when analyzing the stall force of an ensemble of interacting 
molecular motors [28]. To quantify this effect, we therefore define an apparent stall force, 
as the value of force where the velocity drops to less than a small fraction a = 2.5% 
of the value it has for zero force [2S|- In the experimental situation, this bound could 
correspond for instance to the limit of resolution in the velocity measurement. 

The value of the velocity at zero force corresponds to the maximum velocity. When 
F = 0, there is no coupling between the filaments, which behave as independent random 
walkers. The probability to have more than one walker at the leading position is zero 
in the long time limit, which implies q = 1- Therefore, Nq = 1 and the velocity at zero 
force equals the polymerization velocity of a single filament: 

V{F = 0) = d{Uo-Wo), (17) 

which is mainly controlled by the monomer concentration. Now using the expression of 
the velocity at an arbitrary force given by Eq. [TO] the expression of Nq given in Eq. |9] 
and the parametrization of the rates of Eq. [2] for the particular case 7 = 1, we find that 
(^)^fcBT (1 - a){Uo - Wo)N + aUo -{a- q)Wo 

d ^ aUo-ia- q)Wo ' ^ ^ 

Since g <C 1 near stalling, we can write the following more explicit expression 

^S^-^lnfl + --ivl (19) 



'^P^ d ' a y 

In Fig|6| we show the apparent stall force given by Eq. 18 as function of together 



with the theoretical stall force of Eq. 13 
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Nn 



a) 



1 - 



N. 



(20) 



Let us show now that filament condensation at the barrier and the drop in velocity 
occur simultaneously. Assuming for simplicity that 7 = 1, ^ 1 and g ~ in the high 
force regime, we can substitute Eq. [19] into Eq. [9] to obtain: 

Wo 
Uo 

From this we see that since a << 1, the maximum number of filaments at the barrier 
is almost reached. If Vq is the initial velocity and is the finite fraction of filaments 
at the barrier at stall force, we have the equivalence of the following two conditions 

V = aVo^No = {l- a)A'o, 

which shows that filament condensation occurs at the value of the apparent stall force, 
a point which is confirmed by simulations. Indeed, in the case of Fig. |6] the apparent 
stall force is about 12.7 pN, and the condensation visible in Fig. |2]also occurs close to 
12 pN. 

Close to stall force it is also possible to derive an analytic expression for the force- 



velocity relation by substituting into Eq. 10 the expressions of q, given by Eq [14] and 
Eq. [16] Assuming for simplicity 7=1, and using Eq. [11} we obtain with the binomial 
form: 

N-l' 

Uoe-f - Wo 



= A^ 1 



Wo 
Uo 



Wo 
Uo 



(21) 



V = N 1- 



[Uoe-f - Woe 



-N{l-Wo/Uo)-^ 



and with the Poissonian form: 

Wo 
Uo 

When these expressions are expanded close to stall force, one obtains in both cases: 

Wo 



(22) 



5V = N\l 



Uo 



Uoe~f5f. 



(23) 



This indicates an exponential dependence of the velocity close to stalling, which is indeed 
present in the simulations as shown in FigjS] 

To summarize, we have shown in this section that the apparent stall force does not 
scale linearly with A^ as the theoretical stall force but rather as In(A^). The apparent 
stall force is the quantity of experimental interest, it is also near the apparent stall force 
that the condensation transition discussed in a previous section occurs (nothing special 
of that sort occurs near the theoretical stall force). 



4-5. Related experimental work in connection with the model 

In this section we discuss related experimental work. Although a precise comparison 
with the present model is not attempted, we hope that the discussion could be useful 
in identifying some relevant questions in this field. The force generation by parallel 
actin filaments growing out of an acrosome bundle has been measured in Ref . [21] . The 
observation of a plateau in force measurements by optical tweezers is a good indication 
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of the stalling regime, but the measured stall force is very small, comparable with 
that of a single filament, although many filaments are present (about a dozen). These 
results thus stand at odds with the theoretical predictions for the stall force obtained in 
Refs. [151 [20] (and in the present paper). In the present paper, we have emphasized the 
fact that the approach to stalling is slow, which can lead to an underestimation of the 
true stall force. The resolution of the optical tweezers leads to a limit in the detection 
of small velocities, which corresponds roughly to the criterion for the apparent stall 
force used in the previous section. However, with a dozen filaments, the apparent stall 
force should be significantly larger than that of a single filament. Another difficulty is 
that there is no indication in this experiment of the two regimes of low and large forces 
discussed in this paper. At this point, it may be important to say that the results of 
this experiment have not been reproduced, in fact in a new experiment discussed below, 
where the force generated by filaments growing from two magnetic beads outwards, very 
different results have been obtained [22|. In view of all this, we think that the reason 
for these discrepancies may be foundin effects which are not accounted for (such as 
buckling or filament cross-linking) or they maybe attributed more simply to the fact that 
the two experiments have been done in very different biochemical conditions. Indeed, 
the authors of Ref. pi] have used profilin to suppress spontaneous nucleation of actin 
filaments, while profilin was absent in [22J. The use of profilin in Ref. introduced 
complications since profilin also modifies the thermodynamics of the system by binding 
to actin monomers, and possibly interferes with ATP hydrolysis during polymerization. 

The mechanical response of an actin networks confined between two rigid fiat 
surfaces has been probed using a surface force apparatus (SEA) in Ref. [29], and using 
an atomic force microscope (AFM) in [3Dj- Both experiments reported a load history 
dependent mechanical response, which presumably refiects a complex interplay between 
buckling and polymerization forces. This complex interplay makes it difficult to isolate 
the true contribution of polymerization forces. More recently, C. Brangbour et al. 
devised a new experimental setup in which actin is nucleated from magnetic beads which 
are covered by gelsolin [22]. A magnetic field is used to counteract the polymerization 
force, which allows to measure the force- velocity curves. As mentioned above, the results 
of Ref [2T] for the stall force of a single filament are not confirmed: on the contrary, 
the stall force which is obtained is of the order of 40 pN, which corresponds according 
to Eq. 13 to about 25 active filaments. The general shape of these force-velocity curves 
is similar to the ones obtained in this work, but some deviations are present at low 
and high forces. These discrepancies suggest that our model may be too simple to 
fully explain this experiment, and that other aspects may be important. First, it would 
be necessary to go beyond the parallel organization of the filaments in order to better 
model the experimental geometry of Ref. [22] • Secondly, it is probably important to 
account in the model for the possibility of nucleating new filaments from existing ones 
[31] . Thirdly, buckling forces could play an important role in the experiment. Some 
of these effects have been included in previous numerical simulations of branched actin 
networks [I8l[32], but they are typically difficult to study with analytical models of the 
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kind presented here. 
5. Conclusion 

In this paper, we have provided a new theoretical framework to describe the dynamics 
of an ensemble of parallel filaments with no lateral interactions, which are exerting a 
force against a movable barrier. The special cases = 1 and N = 2 can be solved 
exactly, unlike the general case for arbitrary A^, for which we have constructed a 
mean-field approach. We identify two regimes: a non-condensed regime at low force 
in which filaments are spread out spatially, and a condensed regime at high force in 
which filaments accumulate near the barrier. The transition occurs near the apparent 
stall force where the velocity approaches zero. We find that for large N this regime 
where velocity approaches zero occurs at forces significantly lower than the theoretical 
stall force, given by N times the stall force of one filament. In fact, the apparent stall 
force does not scale linearly with N as the theoretical stall force does; instead it scales 
logarithmically. 

On the theory side, several extensions of our work are worth investigating. For 
instance, bundles can be formed experimentally by growing filaments in the presence of 
specific proteins which cross-link the filaments. To describe such a situation, it would be 
necessary to include lateral interactions. Another direction would be to explore the role 
of load sharing, as done in [26] for instance. Although the dynamics will be different, 
we still expect a condensation transition to be present in this case. 

In the end, our model offers a very simplified view of the problem of force generation 
by actin filaments, but precisely for this reason we hope that it can be a useful starting 
point for more refined studies. 

Acknowledgments 

The authors would like to thank A.B. Kolomeisky, P. Sens, R. Pandinhateeri for 
stimulating discussions, and J. Baudry for a careful reading of the manuscript. K. 
Tsekouras would also like to thank J. Elgeti for his help with computational issues. 
This work has been supported by the ANR (french national research agency) under 
contract ANR-09-PIRI-0001-02. 

Appendix: Exact solution of the master equation for case (A^ = 2) 

For this problem, we can derive the following master equation satisfied by p{n, t) the 
probability that there is a gap of n monomers between the two filaments. For n > 2, 



dp{n, t) 
dt 



{U+Wo)[p{ri-l,t)-p{ri,t)] + {W+UoMn+l,t)-p{ri,t)l{K.l) 



and otherwise. 



dp{l.t) 
dt 



2([/+iyo)p(0, t)-{U+W,)p{l, t) + {W+Uo) [p(2, t)l (A.2) 
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dp{0,t) 



dt 



{W + Uo)p{l, t) - 2{U + Wo)p{0, t). (A.3) 



Solving Eq. A.l| at steady state results in the recursion 

{W + Uo)p{n + 1)-{U + W + Uo + WMn) + {U + - 1) = 0, 

which yields two solutions namely 1 and h = ^q^- This means that for n > 2, 
p{n) = p(2)6"-2. 

Using the normalization condition: p{0) + + X]n>2P(2)^"~^ = 1; we obtain 
p(2) = (1 — — p{0) — p{l))- Solving Eq. A.3 at steady state results in: 

U + Wn 

p(l) = 2p(0)^^=2p(0)6. 



Substituting this expression into Eq. |A.2| at steady state yields: 
p{2) = 2p(0)6l 

Equating the two expressions for p{2) and using the expression for p{l) in terms of p{2), 
it follows that 

, , 1-b 
p(0) 



1 + 6 

The probability of having a gap of zero monomers is the probability of having both 
filaments at the barrier; it thus obeys p{0) + g = 1, since q is the probability of having 
only one filament at the barrier. In the end, we find 

+ . (A.4) 

The average number of filaments at the barrier is the sum of q plus twice 1 — q, since 
1 — g is the probability of having both filaments at the barrier. Therefore: 

No = q + 2{l-q) = 2-q, 

and if we substitute q from above we find 

2(W + [/„) 

= u + w + u„ + w,; 
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Figure 1. Representation of the filaments pushing on a barrier (the white vertical 

rectangle on the right, which exerts a force F on the filaments). The right figure 
corresponds to the case that only one filament is in contact with the barrier while 
the left figure corresponds to the case where several filaments are in contact with the 
barrier. The on and off rates of monomers onto free filaments arc Uq and Wq- The 
on-rate on filaments in contact is U, and the off-rate is W when there is only one 
filament in contact and Wq otherwise. 
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Figure 2. Illustration of the condensation transition of actin filaments against a 
barrier. Top: Average barrier velocity vs. force. Symbols represent simulation results, 
the dotted line represents mean- field predictions based on Eq. [TOj Bottom: Average 
number of filaments in contact with the barrier. Symbols represent simulation results, 
the dotted line represents mean-field predictions based on Eq. [9j For both plots, the 
parameters are N — 100, 7 = 1 and C — 0.24/iAf. Inset, Left: Density profile in 
the non-condensed regime (bars) as function of the distance to the barrier, together 
with mean-field theory prediction (line) from Eq. [7j for an applied force F = 2pN, 
which is low with respect to the apparent stall force. Inset, Right: Density profile 
in the condensed regime (bars) as function of the distance to barrier, together with 
mean-field theory prediction (line) from Eq. [7j for an applied force F — 12pN, which 
is close to the apparent stall force of « 12.5pA^. 
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Figure 3. Comparison between tlieoretical and numerical estimates for the parameter 
q, which represents the probability that there is a single filament in contact. Symbols 
represent simulation results, dotted line corresponds to Eq. [16] and continuous line 



corresponds to Eq. 14 (both expressions are mean-field approximations valid in the 



high force regime). The parameters are = 10, 7 = 1 and C = 0.24/iM. 
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Figure 4. Probability distributions of the number of filaments in contact with the 
barrier at various forces. The parameters are N = 100,7 = 1 ^'^d C — 0.24/iM. 
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Figure 5. Average barrier velocity in logaritlimic scale as function of the force in linear 
scale. Note that the velocity decreases to near zero exponentially when approaching 
the theoretical stall force, which is shown with the arrow. For this value of the force, 
the numerical velocity is not strictly zero but it is close to the uncertainty intrinsic to 
the simulation, which is here of the order of lO^^nm/s. The parameters are N = 10, 
7=1 and C = 
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Figure 6. Theoretical stall force F^/^;; (straight line - calculated from Eq. 13 1 and 
apparent stall force, both as computed fr'om simulations Fsimuiations (black symbols) 
and from the mean-field approximation given in Eq. [is] {Fj^pjroximated ~ dotted line) 
vs. number of filaments N . The parameters are 7 = 1 and C — 0.24/iM. 
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Table 1. Parameters characterizing an actin filament barbed end. Wq is the free 
filament depolymcrization rate, ko is the rate constant entering the free filament 
polymerization rate Uq = koC, where C is the concentration of free monomers, d 
is the monomer size and Cc is the critical concentration. 

Wo{s-^) koifiM-^s-^) d(nm) Ce(/xM) 



1.4 



11.6 



2.7 



0.141 



